Dynamics of a one-dimensional spinor Bose liquid: a phenomenological approach 
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The ground state of a spinor Bose liquid is ferromagnetic, while the softest excitation above the 
ground state is the magnon mode. The dispersion relation of the magnon in a one-dimensional 
liquid is periodic in the wavenumber q with the period 2nn, determined by the density n of the 
liquid. Dynamic correlation functions, such as e.g. spin-spin correlation function, exhibit power-law 
singularities at the magnon spectrum, u — > cu m (q, n). Without using any specific model of the inter- 
particle interactions, we relate the corresponding exponents to independently measurable quantities 
duim/dq and dio m /dn. 
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Bosons with an internal degree of freedom, "spin" , ex- 
hibit a ferromagnetic ground stated. The presence of the 
internal states yields an excitation, magnon, in addition 
to conventional waves of the mass density^. Because 
of the ferromagnetic ground-state order, the dispersion 
relation Lu m (q) of magnons at small wave vectors q is 
quadratic. In case of spin-isotropic liquid, this is the 
softest excitation of the system, since ui m (0) = 0, while 
the density waves in the q — ► limit propagate with a 
finite sound velocity v. Away from the isotropic limit 
^m(O) > 0, still magnon is the lowest energy spin exci- 
tation. One may think of a magnon as of a "quantum 
impurity" with the spin opposite to the majority direc- 
tion, moving in a spin-polarized host liquid. 

A spin flip in a compressible system may excite a wave 
of density in it, which affects the spin-spin correlation 
function A m (q,u>). This back-action of the medium is 
the strongest in a one-dimensional (ID) system. In a ID 
Heisenberg ferromagnet on a rigid lattice, A m (q,u>) cx 
/ (q)6 (to — aj m (q)) with f(q) being periodic in the recipro- 
cal lattice. A finite compressibility of liquid transforms 
the J-function response into a power-law edge singularity, 
A m (q,uj) cx (w - u m (q)Y m 6{w - u m (q)), at the magnon 
spectrum. While the existence of that singularity can be 
argued on the basis of the scale invariancei, the exponent 
H m {q) has been evaluated thus far only for an SU(2) sym- 
metric system in the limit of strong repulsion between the 
bosons^, and for an inegrable Yang-Gaudin model^. 

The singularities in A m (q,uS), along with the pre- 
viously studied singularities in the dynamic responses 
of single-component quantum fluids^, and the spec- 
tral function singularity in the quantum impurity prob- 
lem share common root and can be related ^ 10 ' 11 to 
the physics of Fermi edge singularity. The exponent 
Mm (q) > — 1 originates in the interaction of the "quantum 
impurity" with low-energy density waves of the majority 
spin polarization. 

Here we express the exponent /J> m (q) through the dis- 
persion of the magnon mode u> m = u> m (q,n) and its 
derivatives over the wavenumber q and the density of the 
majority spin component n. For this purpose, we find the 
constants of the "quantum impurity" Hamiltonian which 



describes the dynamics of the magnon excitation and its 
interaction with fluctuations of density. That, in turn, 
reduces the problem of finding the critical behavior of a 
dynamic response to a much simpler problem of evalua- 
tion of an excitation energy as a function of n and q. Our 
results for fi m {q) are not limited to small momenta and 
independent of the interaction strength and the presence 
of the SU(2) symmetry. In addition to A m (q, ui), we con- 
sider also the single-particle spectral function Ad{q,uj) of 
a boson with the spin opposite to the majority polariza- 
tion added to the fully polarized liquid. It also displays 
a power-law threshold behavior, and we find the corre- 
sponding wavenumber-dependent exponent Hd(q)- 

To find the exponents fjt m (q) and /J,d(q) we start with 
introducing a representation suitable for the description 
of the low-energy spin dynamics. For definiteness we con- 
sider spin-1/2 particles, and assume all spins point in the 
positive- z direction in the state with maximal spin po- 
larization. Addressing the spin density operator, 



5(x — Xj), 



(1) 



we concentrate first on its z-component. Within the set 
of states of maximal spin polarization, spin density fluc- 
tuations are proportional to the fluctuations of the num- 
ber density. Thus we represent the long-wavelength part 
of density s z (x) in terms of boson field 4> as s z (x) — 
d x <f>/(2ii), where field <j> has the meaning of the displace- 
ment field of spin-up particles. The long- wavelength dy- 
namics of the displacements is described by the Luttingcr 
liquid Hamiltonian^, 
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v 

27 



dx 
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where 9 is the conjugate to <f> variable, [<fi(x), d x 9(y)) — 
iir5(x — y). For Galilean invariant system the Luttingcr 
parameter is K = im/ (m-\v) with m-f being the bare mass 
of spin-up bosons. 

Considering the operator S—(x) in the sub-space of ex- 
citations with energies close to the magnon energy uj m (q) , 
we replace the general form Eq. ([T]) by an operator 



s-(x) ocd^(x)e iqx e i6( - x \ 



(3) 
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Acting on a fully polarized state, this operator "extracts" 
a particle with spin-up, and replaces it (by acting with 
the creation operator S on vacuum) with a spin-down 
particle at the same point x. The representation Eq. ([3]) 
of S-{x) operator is adequate for \q\ < 7m, and requires^ 
a generalization (considered later in this paper) for larger 
momenta. 

The Hamiltonian describing the magnon dynamics 
should include the dynamics of the e?-quasiparticle and 
its interaction with the fields </> and 9 representing the 
density fluctuations in the spin-up state. Similar to 
Refs. MHIIIHj we present the effective "quantum impu- 
rity" Hamiltonian in the form 

H = H + IdxS [u m (q,n) - v m (q)(-id x )]d 



dx[V4,(d x <t>) + V e (d x 6)}dU, (4) 



where v m (<?) = d q uj m (q, n) . The first of the two added 
parts to H establishes ui m (q, n) as the lower edge for the 
spectral function^ 



A m (q,uj) 



]8(u-E v (q))\(v, q \8-(q)\0)\ 



where \u, q) is a many-body eigenstate of the system with 
momentum q and energy E v (q). The gradient expansion 
used in that part is sufficient for finding the behavior of 
A m (q,uj) in the vicinity of the edge. The last term in 
the r.h.s. of Eq. Q describes interactions of the quan- 
tum impurity with the density waves. The strengths of 
such interactions, V<j, and Vg, can be expressed through 
independently measurable characteristics of the system. 

To determine Vg , we note that the corresponding term 
in the Hamiltonian Eq. (j4|) is nothing but a modifi- 
cation of the energy-momentum relation, wj^(g, n) — 
w m ((j) + Vgm^u, for the d-quasiparticle (a mobile "im- 
purity") in the presence of finite velocity u = d x 6/m-\ of 
the fluid formed by the spin-up particles. That allows^ 
one to use the Galilean invariance to find Vg. Indeed, 
in the presence of a flow of the liquid with velocity u, 
the magnon energy remains unchanged in the co-moving 
frame. In the laboratory frame, the magnon momen- 
tum is q + m^u, while its energy is u)^(q + m\u,n) — 
ui m (q,n) + qu + m^u 2 /2. Changing the momentum in 
the last formula, q + m^tt — > q, and comparing the two 
above expressions for u!^(q,n), we find in the limit of 
small u: 



Ve(q) = 



dqU m (q,n). 



(5) 



Hereinafter we assume the mass of bosons is spin- 
independent, m-f = m| = m, thus reducing Eq. |(5j) to 



Vg(q) 



q_ 

rn 



v m (q)- 



(6) 



At \q\ <C 7rn, the magnon spectrum is quadratic in q, i.e. 
oj^q) w w m (0) + q 2 /(2m*) with some effective mass m*; 
in this limit, Vg = [(to* — m)/m]v m . 



To determine V^, we may consider the effect of a long 
wavelength density variation, bp — (l/7r)V0 on the en- 
ergy of the system. According to Eq. (j4]) , it adds a term 
V/ftTrSp to the energy density. This variation should be 
equal to the corresponding value, (duj m /dn) + (dp/dn) 
defined phenomenologically (here p is the chemical po- 
tential of the spin- up bosons). Expressing dp/dn in 
terms of K, one findsi^ 



mi) = — 

7T on 



v 

K' 



(7) 



As follows from the time-reversal symmetry of the prob- 
lem, Vg(q) and V$(q) are correspondingly odd and even 
functions of the momentum q. 

It is instructive to consider separately the q = limit. 
An introduction of a static impurity (spin-down boson 
with q = 0) creates displacement —(ttKV^/v) sgn(y — x) 
in field 0, while operator exp[i#(x)] in the definition of the 
spin density Eq. Q creates a shift 7rsgn(?/ — x). The sum 
of the two is related to the dilatation Si caused by the spin 
flip, -ttKV^/v +tt = w6i, yielding = (v/K)(l - Si). 
Comparing that with Eq. ([7]) at q = 0, we find 



Si = 



K du) m (0,n) 
ttv dn 



Derived from considering a single spin flip, this relation 
under some assumptions can also be derived from the 
thermodynamics of a system at a constant pressure. One 
must assume that magnons do not form bound states, so 
u> m (0, n) = d£/dS with £ being the ground-state energy 
and S the total spin. In the presence of SU(2) invariance 
d£/dS = 0, thus flipping a spin with q = does not 
cause dilatation, Si = 0. 

Once the interaction constants Vg and V^ are estab- 
lished, we may proceed in full analogy with Ref. [l3l |: 
(i) re-scale variables, = 4>Vk, 6 = Bj\[K thus giving 
Eq. ([2]) the appearance of a Hamiltonian of free "par- 
ticles"; (ii) eliminate the linear in cj), 8 part of the full 
Hamiltonian, Eqs. ((2]) and (j4|), by a unitary transforma- 
tion with proper values of S+, 5—, 



-i/<fa{^[9"W-#W]- 



-[6i(x)+^(a)]}d(a;)d t (x) 



(8) 

(iii) express the exponents of the sought correlation func- 
tions in terms of the phase shifts S+ and 5- which "quan- 
tum impurity" (d)d = 1) causes for co-moving (+) and 
counter-propagating (—) "particles". The values of the 
phase shifts can be written in a relatively compact form 
as: 



1 



v m (q) T v 



0) 



We are now in the position to evaluate the correlation 
functions of interest. Using Eq. §3§ we may represent 
A m (x, t) as 

A m (q,u) (10) 
r dxdte lqx - luJt (0\d(0,0)e-' wio ' 0} e l9(x ' t) d^x 1 t)\0). 
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In the case of spinor Bose liquid one may envision tun- 
neling of a spin- up or a spin-down particle into otherwise 
fully spin-up polarized system. The spectral function for 
a spin- up particle is identical to the one evaluated for the 
liquid of spinless bosons^; the corresponding tunneling 
theshold spectrum is uj = vq. For a spin-down parti- 
cle, the tunneling threshold is determined by the magnon 
spectrum, and the tunneling probability is proportional 
to the spectral function of the <i-quasiparticle, 



Ad{q,u) 



dxdt i 



Jt (0\d(x, t)tft(o,o)|0) ■ (11) 



Functions A m (q,u>) and Ad(q,ui) exhibit power-law be- 
havior above the threshold, 



A m . d (q,uj) oc Q[u - oj m (q)] [w - w m (q)] 



A»m,d(g) 



Using the transformation Eq. JHJ) and the standard meth- 
ods of bosonisation, one finds for fid(o) in the region 
\q\ < nn: 



2- 



2tt 



like in Ref. [13|. Similarly, the exponent /U m (g) at \q\ < 
nn reads as 



1 



2y/K ^ 2tt 



1 



*-(g) 
2y/K 2tt 



(13) 

Equations QTZty and (13]) together with Eq. © relate the 
exponents of the correltaion functions to the properties 
of the magnon branch of excitation spectrum u> m (q,n) 
for the principal interval of momenta \q\ < nn. 

Due to the peculiarity of ID systems, the lowest-energy 
excitations corresponding to a single flipped spin or to 
an added spin-down particle at a given momentum q, 
are periodic functions of the momentum, Lu m (q,n) — 
uj(q — 2nnl,n) for any integer I. To extend the above 
results for fi m ,d(q) beyond the principal interval of mo- 
menta, we notice that introducing a d-quasiparticle with 
the lowest energy amounts to a momentum boost 2nnl 
of the spin-up liquid and exciting a magnon with the 
residual momentum belonging to the principal interval 
—nn < q — 2nnl < nn. The boost accompanying the 
creation of a d-quasiparticle corresponds to the modified 
definition of the spin-down particle creation operator: 

<ft(x) e- 2il ^d^(x) ; d(x) — d(x)e 2lt ^ . (14) 

Performing these replacements in Eq. (jlip and repeating 
the steps which have lead to Eq. (Tj"2"|) , we find now: 



fJ-d(q) = -1 + 



Mg* 

2tt 



Mg* 

2tt 



where q* = q — 2nnl with the integer / chosen to have 
\q*\ < nn. A similar procedure for A m (q, uS) yields 



Mm(g) 



-1 



1 



Mg*) 

2y/K ^ 2tt 

1 Mg*) 
2n 



2VK 



iVk) . (16) 



Equations (Ti"5|) , (TTB]) along with Eq. © provide the values 
of the edge exponents for an arbitrary momentum q. 

Transferring momentum to the liquid as a whole allows 
tunneling of a particle at low energy, oj — > u> m (q* , n) , even 
at high momentum |g| > nn. The price for that, however, 
is a reduced tunneling probability reflected by the pres- 
ence of integer I in Eq. (fT5|l : while the spectrum uj m (q) 
is periodic, the exponent fJ-d{q) is increasing with moving 
from one period to another, with larger |/|. The sup- 
pressed tunneling probability is a manifestation of the 
orthogonality catastrophe. Similarly the exponent /j. m , 
describing the probability of the spin-flip photon absorp- 
tion near the edge, is increased due to the orthogonality 
catastrophe. Indeed, for \l\ > the final state includes 
the spin-flipped particle along with the moving fluid with 
the momentum 2nnl, which has the progressively smaller 
overlap with the initial state of the fluid at rest. 

The periodic dispersion relation reaches its maxima at 
q = nn(2l — 1). Depending on the microscopic interac- 
tion strength between the bosons, the magnon velocity 
v m = d q u) m {q,n) may have jumps at q = nn(2l — 1), or 
be a continuous function. In the latter case, obviously, 
v m (nn(2l—l)) = 0. The transition between the two types 
of behavior upon the increase of the interaction strength 
is equivalent to the "quantum phase transition" in the 
Kondo problem controlled by tuning the exchange con- 
stant through zeroii. The v m (nn(2l — 1)) = regime 
corresponds to the strong-coupling side of the transition. 
The developed above Luttinger liquid representation is 
applicable on either side of the transition, similar to the 
scatterling phase description of the low-energy physics 
of the Kondo problem. The region of applicability in 
lo — u) m (q,n), of course, gets narrow close to the tran- 
sition point, as the corresponding Kondo energy scale 
becomes small. In the strong-coupling regime, v m and 
dn^m{q-,n) have no discontinuities at q — nn(2l — 1). It 
is not clear apriori that the same is true for the expo- 
nents Hm,d(q)' after all, the definition of the response 
functions involves different operators, see Eq. (Ti"4")) . at 
subsequent intervals of momenta. It is quite striking 
to see directly from Eqs. (fT5|) . (fTB"|) and (J9]) the con- 
tinuity of Hd{q) and fJ, m (q). Indeed, substitution of 
y m(g* = inn) = 0, the use of relations K — nn/(mv) 
and d n u m (q* ,n)\ q . =ml = d n u m (q* , n)| g *=_ 7r „ in Eq. © 
yields: 



5±{nn) 8±{~nn) 



2n 



2n 



- Vk. 



(17) 



(15) 



With the help of Eqs. (Tr5|) and (TT6"]) this immediately 
implies that for q — > nn(2l — 1) ± both fj.d(q) and /J, m (q) 
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are continuous functions and their first derivatives over 
q are continuous as well. At the "weak-coupling" side of 
the transition, the exponents, together with v m (q), are 
discontinuous at q — im(2l — 1). 

Around the local minima, q — 2ml, i.e. at q* = 0, one 
finds 



where a = — (n/2m*)dm* /dn. Similarly, 



/id(27mZ) 



-1 



(1 



and 



H m {2iml) = — 1 



2K 



2K 



2l 2 K 



2VK. 



(18) 



(19) 



The fact that Hd{q = 0) 7^ — 1 is due to the orthogonality 
catastrophe^: the tunnelled spin-down boson shakes up 
the liquid of spin- up particles. The same mechanism have 
led 8 to non-trivial tunneling exponents fj,± for spinlcss 
bosons in the Lieb-Liniger— model, once the exponents 
are evaluated beyond the Luttinger liquid approximation. 
The direct comparison of Eq. (|18|> with the corresponding 
result, Eq. (22) of Ref. [1] is possible only at if — > 1, 
when the shake-up becomes independent of the impurity 
velocity relative to ±u; the two equations agree with each 
other, yielding fi^ — ~fl± = 1/2. The similar physics is 
at work for A m (q, u>) correlation function away from the 
SU(2) symmetric point, i.e. when Si =/= 0. The amplitude 
of a spin-flip process for a finite dilatation Si is suppressed 
by the orthogonality, resulting in /i m (0) 7^ — 1. 

In the SU(2) invariant case, where duj m (0,n)/dn cx 
Si = 0, one may also deduce universal results for the 
momentum dependence of /i m ,d(q) in the region of small 
momenta \q\ <C m*V, 



1 Kq 2 
2K 2{irn) z 



1 + (2 + 2cr) , (20) 

to 



Mm(e) = -1 



Kg 2 
2(im) 



4(7 ■ 



The q 2 term of the /i m (<?) exponent depends only on the 
parameters of the Luttinger liquid and for the strong cou- 
pling limit was derived in Refs. It is interesting that 
the q 4 term here may be expressed through the effective 
mass of the magnon mode to* . The latter was evaluated 
in various limits for the integrable contact-interaction 
models, which allows us finding a in Eqs. (|2"0"]) and 
(f2"Tj) . In the strong coupling limit to* — 37to/(27t 2 ), 
where 7 = mg/n 1 and thus a = 1/2; here g is 
the interaction strength. In the weak coupling limit 
to* = m(l + 2^/7/371-) and thus a = ^f/Sir < 1. These 
considerations are also applicable for the q 2 term in fid(q)- 

In conclusion, the dynamic response functions 
Eqs. (fTTj|) . (fTl"]) of a homogeneous ferromagnetic one- 
dimensional Bose liquid exhibit power-law asymptotes at 
the threshold defined by the spectrum of the magnon 
uj m (q). Independent of any model assumptions, the cor- 
responding exponents n m (q) and l^d(q) at any wave vec- 
tor can be expressed in terms of a few independently 
measurable parameters: the sound velocity v, the corre- 
sponding Luttinger liquid parameter K, and the deriva- 
tives of u m with respect to q and liquid density n. Fur- 
ther simplification of /J. m (q) and fJ,d{q) is possible in the 
vicinities q = 2irnl of the minima of magnon spectrum. 
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